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ABSTRACT 

Using XMM we detect faint unresolved X-ray emission from the Compact 
Steep Spectrum radio galaxy 3C 303.1. We detect a thermal component at 
kT ~ 0.8 keV which seems likely to be produced in the ISM of the host galaxy. 
There is evidence for a second component in the spectrum whose nature is cur- 
rently ambiguous. Plausible hypotheses for the second component include (1) 
hot gas shocked by the expansion of the radio source, and (2) Synchrotron self- 
Compton emission from the southern radio lobe if the magnetic field is below the 
equipartition value by a factor of ~ 3.5. 

Subject headings: galaxies: active — galaxies: jets — galaxies: individual 3C 303.1 
— X-rays: galaxies 

1. INTRODUCTION 

The progenitors of the large-scale powerful classical double (3CR FRII) sources seem 
likely to be the GHz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) radio 
sources (e.g., O'Dea, Baum, & Stanghellini 1991; Fanti et al. 1990; Fanti et al. 1995; Read- 
head et al. 1996a,b; O'Dea & Baum 1997, for a review see O'Dea 1998). The GPS and CSS 
sources have generally simple and convex radio spectra with peaks near 1 GHz and 100 MHz, 
respectively. 

Models for the evolution of powerful radio galaxies (Fanti et al. 1995; Begelman 1996; 
Readhead et al. 1996b; O'Dea & Baum 1997; De Young 1997; Kaiser & Alexander 1997; 
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Kaiser, Dennett-Thorpe, & Alexander 1997; O'Dea & Baum 1997; Blundell, Rawlings & 
Willott 1999; and cf Snellen et al. 2000) are generally consistent with a scenario in which 
these sources propagate from the ~ 10 pc to Mpc scales at roughly constant velocity through 
an ambient medium which declines in density while the sources decline in radio luminosity. 

There is increasing evidence that the GPS and CSS sources interact with the host galaxy 
as they propagate through it. Broad- and narrow-band HST imaging of CSS radio sources 
(De Vries et al. 1997, 1999; Axon et al. 2000) have shown that CSS radio galaxies at all 
redshifts exhibit emission-line gas which is strongly aligned with the radio source. The 
alignment is much stronger than seen in low redshift large-scale radio galaxies and is similar 
to that seen in high redshift (z j> 0.6) radio galaxies. The close association between the gas 
and the radio source suggests that the latter interacts strongly with the ambient medium as 
it propagates through the ISM. 

In 3C 303.1 cooling time arguments (e.g., Pedlar, Dyson, & Unger 1985; Taylor, Dyson, 
& Axon 1992) applied to the gap between the radio hot spots and the onset of bright emission- 
line gas are consistent with lobe expansion velocities J> 6000 km s _1 for ambient densities 
of n ~ 1 cur 3 (De Vries et al. 1999). The broad and highly structured spatially integrated 
[OIII]A5007 line widths observed by Gelderman & Whittle (1994) strongly suggest that 
the radio source is dominating the emission-line kinematics. HST/STIS long-slit medium- 
dispersion spectroscopy of the aligned emission-line gas has shown that the kinematics of 
the gas are consistent with cloud motions driven by interaction with the expanding radio 
source (O'Dea et al. 2002). Constraints on the bow shock velocity from lifetime and cooling 
arguments are in the range few x 10 3 km s _1 to few x 10 4 km s _1 (O'Dea et al. 2002). Analysis 
of the low dispersion STIS spectra reveals that the ionization diagnostics are consistent with 
strong contributions from shocked gas (Labiano et al. 2005). HST UV imaging shows a clear 
detection of UV continuum light tightly aligned with the radio source (Labiano et al. 2006). 
These observations are all consistent with the radio source in 3C 303.1 interacting with gas 
clouds in its environment, shocking them and triggering star formation. 

The interaction of the radio source with its environment may also produce observable 
signatures in the X-ray. Heinz, Reynolds & Begelman (1998) suggested that compact evolv- 
ing radio galaxies might create shells of hot shocked gas as they expand into the ambient 
medium. X-ray observations have found examples of both hotter (e.g., Kraft et al. 2003; 
Wilson, Smith & Young 2006) and cooler (e.g., Nulsen et al. 2005; Soker, Blanton & Sarazin 
2002; Blanton et al. 2001) shells of ambient gas swept up by expanding radio galaxies. 

In this paper we present XMM-Newton X-ray observations of the CSS radio galaxy 
3C 303.1. We discuss the results and their implications for our understanding of CSS radio 
sources. 
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2. THE XMM OBSERVATIONS AND RESULTS 

The properties of the source are given in Table 1. We observed 3C 303.1 with XMM- 
Newton for about 40 ks on August, 18, 2003. Three EPIC cameras, MOS-1, MOS-2 and pn, 
were operated with medium optical blocking filters in "Prime Full Window" mode, which 
covers the full field-of-view (FOV) of 30' diameter. Because a fraction of the flux collected by 
the telescopes hosting the MOS cameras is directed to dispersing grating arrays, the photons 
collected by each MOS camera are fewer than those collected by the pn camera. 

The data were reprocessed using the XMM-Newton Science Analysis System (SAS) 
v6.1.0. EMCHAIN and EPCHAIN were employed to obtain the photon event lists. Then we 
used EVSELECT to select single, double, triple and quadruple (PATTERN 12) for MOS, 
and single and double events (PATTERN ^ 4) for pn with the energy range from 0.2 keV 
to 12 keV. FLAG was set to zero to reject the events which are close to the CCD gap and 
bad pixels. Because of the high EPIC background due to contamination by solar protons we 
restricted the data using Good Time Interval (GTI) filters. Details of the resultant exposure 
times and count rates in the three cameras are given in Table 2. The pn image, showing 
the detection of a faint source at the location of 3C 303.1 is given in the left-hand panel of 
Figure 1. 

2.1. Source Size 

The radial profiles of the source counts are consistent with the X-ray source being 
unresolved by the pn and MOS detectors. This is consistent with the emission being produced 
by the radio source (1.7 arcsec) and/or the host galaxy (20 kpc = 4.9 arcsec). Galaxy counts 
by Harvanek et al. (2001) show no evidence for a cluster of galaxies surrounding 3C 303.1. 

2.2. Spectral Fitting 

We chose an on-source extraction circle of 40 arc-second radius, which at 1.5 keV corre- 
sponds to ~ 87% encircled energy for the MOS and ~ 83% encircled energy for the pn. The 
MOS local background region was an annulus surrounding 3C 303.1 with inner and outer 
radii of 40 and 80 arcsec, respectively, with two unrelated sources masked. For the pn the 
background region was a circle the same size as the on-source region but from a source-free 
region of the same RAWY. For each camera the task ESPECGET was employed to generate 
a background-subtracted spectrum, grouped with 25 counts minimum per spectral bin to 
improve the statistics, a Redistribution Matrix File (RMF), and an Auxiliary Response File 
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(ARF). Data above 8 keV and below 0.4 keV were excluded from the spectral model fitting 
due to the poor statistics at these energies. The resulting count rates are given in Table 2. 

We fitted the pn and MOS spectra simultaneously with several models consisting of dif- 
ferent combinations of absorption with power law (PL) and Raymond-Smith (RS) thermal 
models using XSPEC 12.2.0cf. The results are given in Tables 3-5. We obtain an accept- 
able fit for a RS model with Galactic absorption (frozen), solar abundance (frozen), and a 
temperature of kT~ 0.8 keV (Table 3). We do not obtain an acceptable fit with a single PL 
model with either Galactic absorption alone or excess intrinsic absorption. Our single RS 
fit improves if we add a PL component (Table 4), though the photon index is very poorly 
constrained by our data. Assuming a plausible photon index, e.g., T = 1.5, there is no need 
for absorption in excess of Galactic towards the PL component. In order to constrain a 
possible hot component due to gas shocked by the radio source (§3.2) we considered a two 
temperature model with Ti = 0.8 and T 2 = 45 keV (Table 5). Such a model is also an 
improvement over a single temperature thermal model. Figure 2 shows the spectral data 
together with the model components of Table 4 folded through the detector response files. 

3. DISCUSSION 

3.1. The Thermal Component 

The luminosity (L x ~ 2 x 10 42 ergs s _1 ) and temperature (kT~ 0.8 keV) of the RS 
component are consistent with the values found in elliptical galaxies hosting radio sources 
(e.g., Worrall & Birkinshaw 2000) and possibly small groups (e.g., Helsdon & Ponman 2000). 
However, the X-ray emission is both too faint and too cool to be consistent with a cluster 
of galaxies (e.g., Arnaud & Evrard 1999). This suggests that we have detected the ISM of 
the host galaxy or possibly the intragroup medium (if there is a galaxy group). 

3.2. Constraints on shock heated gas 

The expanding radio source will run over and shock the ambient gas possibly producing 
detectable emission (e.g., Heinz et al. 1998). We calculate the expected properties of the 
shocked gas following Worrall et al. (2005). The sound speed in the 0.8 keV ISM is c s ~ 
460 km s -1 . If the propagation velocity of the radio source is as high as 6000 km s _1 as 
suggested by the cooling-time arguments for the emission line nebula (De Vries et al. 1999), 
the Mach number is M ~ 13. For such a strong shock the Rankine-Hugoniot conditions 
give a temperature ratio between shocked and pre-shocked gas of T 2 /Ti ~ 54, which implies 
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T 2 ~ 43 keV. We note that a two-temperature model with T\ = 0.8 and T 2 = 45 keV is 
permitted by the data (Table 5). 

Taking into account the uncertainties, Table 5 shows that the normalizations of the 
heated (shocked) to cool (unshocked) gas is a factor of roughly 0.9 to 6. The 1-keV 
normalization is proportional to n 2 V, where n is the number density of electrons or pro- 
tons and V is the volume. Since the Rankine-Hugoniot conditions for a strong shock give 
"-shocked /Vnshockcd = 4, the X-ray spectral results are consistent with V mishockcd / V shockcd ~ 3 to 
18. What is the volume of the shocked gas? We detect radio emission from the cocoon but 
not from the hot shocked gas which lies between the cocoon and the bow shock. We use the 
numerical simulations of Carvalho & O'Dea (2002) as a guide to estimate the ratio of the vol- 
umes of hot shocked gas and cocoon. The estimated ratio is model dependent and varies with 
the jet density contrast and Mach number. For a light, supersonic jet the numerical simula- 
tions suggest a value of about 3 for the ratio of hot shocked gas to cocoon volume. Adopting 
a value of 3 along with the estimates of cocoon volume from the MERLIN observations gives 
a volume of shocked gas of about 14 kpc 3 . We find the volume of unshocked gas is large 
enough to be consistent with a realistic galaxy atmosphere if Kinshocked/ ^shocked is towards 
the upper value permitted by the X-ray spectral fitting. For example, if Kmshocked/ ^shocked 
is 18, a spherical galaxy of uniform density would be required to have a radius of about 4 
kpc. Such a galaxy has a similar n 2 V to one of the same central density modeled with a (3 
model of (5 = 2/3 and core radius of ~ 3 kpc. Chandra measurements of the atmospheres 
around relatively powerful nearby radio galaxies are fitted to (5 models that typically bracket 
these parameter values (e.g., Kraft et al. 2005, Worrall & Birkinshaw 2005). High- resolution 
X-ray imaging spectroscopy would probe the density and temperature structure and should 
confirm whether such a high-temperature component exists. 



3.3. The Power-Law Component 

Here we consider the possibility that there is a significant non-thermal power-law com- 
ponent, i.e., the second component in the fit to the spectrum is not due to a hot shocked 
gas. The lack of a high intrinsic absorbing column towards the PL component suggests 
that it is not from accretion-disk emission that would be expected to be obscured in this 
powerful two-sided radio galaxy, and the fact that the nucleus is undetected both at 1.6 
GHz (Sanghera et al. 1995) and at 5.0 GHz (Liidke et al. 1998) disfavors the X-ray emission 
originating from a small-scale jet. X-ray binaries produce emission that can mimic a hard 
power law, but the luminosity of > 10 42 ergs s _1 is one or two orders of magnitude more 
than is typically measured in elliptical galaxies (e.g., Kim & Fabbiano 2004). 
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The most likely explanation for a power-law component would appear to be inverse 
Compton emission from the lobes of the radio galaxy. We estimate the expected flux (e.g., 
Harris & Grindlay 1979; Harris & Krawczynski 2002; Hardcastle, Birkinshaw & Worrall 
1998) using the MERLIN radio measurements of Sanghera et al. (1995), where at 1.6 GHz 
the Southern lobe contains about 90% of the total radio flux density. We find that, because 
of the small size of the lobes, Synchrotron self-Compton (SSC) dominates the X-ray yield 
from scattering of the cosmic microwave background radiation. 

Results of the SSC calculations are presented in Table 6. In these calculations we adopt 
a single- component power-law electron spectrum of slope consistent with the observed radio 
synchrotron spectrum (Sanghera et al.). In order to maximize the SSC output under the 
equipartition assumptions, we set the minimum electron energy to no lower than is required to 
produce the radio emission, and we assume no non-radiating particles (protons) are present. 
If the magnetic field is at the equipartition strength, the predicted SSC X-ray flux is too 
low by factors of ~ 350 in the Northern lobe and ~ 16 in the Southern lobe. However, 
if the field is below the equipartition level by factors of ~ 16 and ~ 3.5, in the N and S 
lobes, respectively, the predicted flux is comparable to that in the measured PL component. 
There is some evidence that the magnetic fields in some radio lobes are indeed about 1/3 of 
the equipartition strength (e.g., Carilli, Perley & Harris 1994; Wellman, Daly & Wan 1997; 
Croston et al. 2005), and thus there is additional justification for the PL component arising 
from SSC emission primarily in the S lobe of 3C 303.1, with a magnetic field about 1/3 of 
the equipartition value. If this is the correct interpretation it suggests that the mechanism 
for keeping the average B field below equipartition occurs in young radio galaxies and is not 
just a feature of large FRIIs. 



3.4. Comparison with other Radio Galaxies 

There are now X-ray observations of a dozen GPS radio galaxies (O'Dea et al. 1996, 2000; 
Iwasawa et al. 1999; Guainazzi et al. 2004, 2006; Vink et al. 2006). 4C+55.16 appears to be 
in a cluster with a cooling flow (Iwasawa et al. 1999). The others are dominated by a power- 
law component with intrinsic absorption of N(H) ~ 10 22 to ~ 10 24 cm~ 2 (Guainazzi et al. 
2006; Vink et al. 2006). This intrinsic absorption is larger than the HI absorption along the 
line of sight to the radio emission in the small number of cases where the radio measurements 
have been made (Vink et al.), thus suggesting that the absorbed X-ray emission arises from 
the vicinity of an accretion disk. 

Large narrow-emission-line FRII radio galaxies generally have two nuclear X-ray com- 
ponents (Belsole et al. 2006). One has intrinsic absorption similar to that measured for 
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most GPS sources, and is interpreted as accretion related, and the other is unabsorbed and 
identified as arising from jet emission within the radio core. Our X-ray measurements of the 
CSS radio galaxy 3C 303.1 find no evidence for absorbed emission - it is either absent or 
very weak. If absent, it may suggest that the central engine of 3C 303.1 is more similar to 
that of an FRI radio galaxy, where absorbed X-ray components are rarely seen (Evans et al. 
2006 - but note the cases of Cen A [Evans et al. 2004] and NGC 4261 [Zezas et al. 2005]). 
The radio core of 3C 303.1 is undetected, and if the second X-ray spectral component is 
indeed non-thermal, then it is much more plausibly associated with the radio lobes than 
with the radio core. The alternative hypothesis is that we are seeing shocked gas. In either 
case this hard emission should be spatially resolved in Chandra observations. It appears 
that the weak core of 3C 303.1 is an advantage in allowing other X-ray components to be 
seen more clearly. 

4. SUMMARY 

Using XMM we detect faint unresolved X-ray emission from the CSS radio galaxy 
3C 303.1. We detect a thermal component at kT ~ 0.8 keV which seems likely to be 
produced in the ISM of the host galaxy. There is evidence for a second component in the 
spectrum whose nature is currently ambiguous. It does not appear to be related to the core, 
giving 3C 303.1 the advantage over other GPS/CSS radio galaxies that extended components 
show up more clearly. Plausible hypotheses include (1) hot gas shocked by the expansion of 
the radio source, and (2) SSC emission from the southern radio lobe if the magnetic field 
is below the equipartition value by a factor of about 3.5. In the former case this would be 
the first detection of strongly shocked gas around an FRII radio source. In the latter case 
it would be a rare instance of the detection of SSC emission from a radio lobe, as distinct 
from a small localized hot spot. Deep Chandra observations could confirm the presence of 
this extended component and improve our understanding of its origin. 
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Table 1. Source Properties 



Parameter 3C 303.1 



ID G 

redshift 0.267 

RA (J2000) 14:43:14.8 

DEC 77:07:28 

scale (kpc/arcsec) 4.07 

radio power lognPsGHz (Watts Hz -1 ) 26.53 

radio source total angular size 6 (arcsec) 1.7 

linear size D (kpc) 6.9 

integrated emission line flux F(OIIIA5007) (10~ 15 ergs s _1 cm~ 2 ) 28 

spectral age (yr) 1 x 10 5 

advance speed (v/c) 0.07 



Note. - We adopt the standard cosmology with H a = 71 km s 1 Mpc 



Q m = 0.3, and = 0.7. The integrated emission line flux is from Gelderman & 
Whittle (1994). The spectral age is estimated by fitting a continuous injection 
model to the integrated radio spectrum and is taken from Murgia et al. (1999). 
The advance speed is estimated using 2 v = linear size / spectral age. 

Table 2. XMM Observations of 3C 303.1 



Parameter pn Mosl Mos2 



Net Exposure Time (ksec) 25.93 32.13 32.67 

Net Count Rate (counts/s) 8.15E-3± 1.92E-3 2.15E-3± 7.22E-4 1.62E-3± 6.97E-4 
Total Counts 628 218 191 
Total Counts 0.4-8 keV 483 127 128 
Net Counts" 0.2-12 keV 211 69 53 



Note. — These numbers are after GTI filtering. Observation date is August, 18, 2003. 
a background subtracted. 



Table 3. Spectral Fits: Single Thermal Component 



Parameter Value 



External Absorption n H (xlO 22 cur 2 ) 0.031" 

Abundance 1.0 fe 

Temperature kT (keV) 0.82 ± 0.05 

Luminosity 0.4-8.0 keV (10 42 ergs s" 1 ) 2.4 
1 keV normalization (10~ 6 photons cm~ 2 sec -1 keV -1 ) 7.7 ± 2.1 

Chi-squared 24.63 

degrees of freedom 27 



Note. - - A fit to wabs*raymond using the energy range 0.4 to 8.0 
keV over 29 bins. "Absorption is frozen at the foreground Galactic 
value. 6 Abundance is frozen at the solar value. 
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Table 4. Spectral Fits: Thermal Component and Powerlaw 



Parameter 




Value 


External Absorption n H (xlO 22 cm -2 ) 




0.031 a 


Abundance 




1.0 6 


Temperature kT (keV) 




0.77 ±0.09 


Luminosity 0.4-8.0 keV (10 42 ergs s _1 ) 




2.0 


1 keV normalization (10~ 6 photons cm 


" 2 sec" 1 keV" 1 ) 


5.9 ±4.3 


Absorption towards the PL % (xlO 22 


cm -2 ) 


0.47 ±1.38 


Photon index T 




1.5 C 


1 keV PL normalization (10~ 6 photons 


cm -2 sec -1 keV -1 ) 


i qo.i 


Luminosity 0.4-8.0 keV (10 42 ergs s" 1 ) 




2.2 


Chi-squared 




17.17 


degrees of freedom 




25 



Note. — A fit to wabs* (raymond ± zwabs*pow) using the energy range 
0.4 to 8.0 keV over 29 bins. "Absorption is frozen at the foreground Galactic 
value. 6 Abundance is frozen at the solar value. c The power-law photon 
index is poorly constrained by our data, so we froze a value of T = 1.5. A 
value of T = 2.0 produces similar results. 
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Table 5. Spectral Fits: Warm and Hot Thermal Components 



Parameter 




Value 


External Absorption nn (xlO 22 cm 




0.031 a 


Abundance 




1.0 6 


Temperature kT (keV) 




0.8 C 


Luminosity 0.4-8.0 keV (10 42 ergs s 




1.7 


1 keV normalization (10~ 6 photons 


cm -2 sec -1 keV -1 ) 


5.3 ±1.7 


Temperature kT (keV) 




45 c 


1 keV normalization (10~ 6 photons 


cm -2 sec -1 keV -1 ) 


13.6 ±7.6 


Luminosity 0.4-8.0 keV (10 42 ergs s' 




2.8 


Chi-squared 




18.77 


degrees of freedom 




27 



Note. - A fit to wabs* (raymond ± raymond) using the energy 
range 0.4 to 8.0 keV over 29 bins. "Absorption is frozen at the 
foreground Galactic value. 6 Abundance is frozen at the solar value. 
c Temperature frozen at 0.8 and 45 keV. 
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Table 6. Estimated Synchrotron Self-Compton Emission 



Parameter 


N Lobe 


S Lobe 


Radio size (mas) 


440 x 200 


180 x 80 


1.6 GHz radio flux density (mJy) 


140 


1435 


radio spectral index a r 


1.1 


1.2 


minimum electron Lorentz factor 7 min 


1000 


400 


Equipartition B field (milliGauss) 


0.17 


0.7 


Predicted SSC X-ray 1 keV flux density (Jy) 


3.6 x 10~ 12 


8.0 x 10" 11 


B field correction factor 


0.06 


0.29 



Note. - The radio data are from 1.6 GHz MERLIN observations pre- 
sented by Sanghera et al. (1995). The B field correction factor is the factor 
needed to reduce the B field so that the predicted SSC emission is consistent 
with the observed PL component. 
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Fig. 1.— 3C 303.1. (Left). XMM pn image showing detection of 3C303.1. A circle of 15 
arcsec radius is shown around the source. (Right). Optical image from the STScI digitized 
sky survey. A circle of 10 arcsec radius is shown around the source. 



,-/\9>l ,-S siunoo paz||BLUJou sienpjssj 



Fig. 2. — The XMM spectrum of 3C303.1 showing the pn (darker) and MOS (lighter) data 
and the model of Table 4. 



